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Abstract 



A factorization theorem for P-wave quarkonium production, recently derived by Bod- 
win, Braaten, Yuan and Lepage, is applied to T ^ Xcj + X, where Xcj labels the ^Pj 
charmonium states. The widths for Xcj production through color-singlet P-wave and 
color-octet 5- wave cc subprocesses are computed each to leading order in a^. Experi- 
mental data on T ^ J /ip + X is used to obtain an upper bound on a nonperturbative 
parameter (related to the probability for color-octet S-wave cc hadronization into P- 
wave charmonium) that enters into the factorization theorem. The bound obtained 
here adds to the limited information so far available on the color-octet mechanism for 
P-wave quarkonium production. 



1 



Factorization theorems play a basic role in perturbative QCD calculations of many 
hadronic processes. A well known factorization theorem for the decay and produc- 
tion of S-wave quarkonium follows from a nonrelativistic description of heavy quark- 
antiquark {QQ) binding Nonperturbative effects are factored into Rs{0), the 
nonrelativistic wave function at the origin, leaving a hard QQ subprocess matrix ele- 
ment that can be calculated in perturbation theory. This factorization is valid to all 
orders in the strong coupling a^, and to leading order in f ^, where v is the typical 
center-of-mass velocity of the heavy quarks. 

Remarkably, the correct factorization theorems for the decay and production 
[0] of P-wave quarkonium have only recently been derived. These new theorems re- 
solve a long standing problem regarding infrared divergences which appear in some 
cases to leading order in the rates for P-wave QQ states In previous phenomeno- 
logical calculations, the divergence was replaced by a logarithm of a soft binding 
scale, such as the binding energy or confinement radius @, However, a rigorous 
calculation requires that one consider additional components of the Fock space for 
P-wave quarkonium, such as \QQg), where the QQ pair is in a color-octet S'-wave 
state, and g is a soft gluon |^, ^. 

A renewed study of the decay and production of P-wave quarkonium is therefore 
of considerable interest, since one may gain new information on a nonperturbative 
sector of QCD that has largely been neglected in the quark model description of 
heavy quarkonium. This is also of practical consequence; for example, J/ip production 
provides a clean experimental signature for many important processes, and P-wave 
charmonium states have appreciable branching fractions to J/ip- 

In this paper the factorization theorem for P-wave quarkonium production is 
applied to T ^ Xcj + X, where Xcj labels the ^Pj charmonium states. The widths for 
Xcj production through color-singlet P-wave and color-octet S'-wave cc subprocesses 
are computed each to leading order in a^. Experimental data on T ^ J/ip+X is used 
to obtain an upper bound on a nonperturbative parameter (related to the probability 
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for color-octet S'-wave cc hadronization into P-wave charmonium) that enters into 
the factorization theorem. The bound obtained here adds to the hmited information 
so far available on the color-octet mechanism for P-wave quarkonium production. 
The color-octet component in P-wave decay was estimated in Ref. from measured 
decay rates of the Xci and Xc2- A rough estimate of the color-octet component in P- 
wave charmonium production was obtained in Ref. from data on B meson decays; 
however, an accurate determination in that case requires a calculation of next-to- 
leading order QCD corrections to the color-singlet component oi B Xcj + X, which 
is so far unavailable 0. 

The factorization theorem for P-wave quarkonium production has two terms, and 
in the case of T decay takes the form: 

r{T ^Xcj + X) = H,r^{T ^ccCPj)+X;fi) 

+ (2 J + l)H',ifi)ts{T ^ ccCSi) + X). (1) 

Fi and Fs are hard subprocess rates for the production of a cc pair in color-singlet 
P-wave and color-octet S'-wave states respectively. The quarks are taken to have van- 
ishing relative momentum. The nonperturbative parameters Hi and Hg are propor- 
tional to the probabilities for these cc configurations to hadronize into a color-singlet 
P-wave bound state. Hi, H^ and Fg are independent of the total angular momentum 
J. 

This factorization theorem is valid to all orders in and to leading order in 
f^. The hard subprocess rates are free of infrared divergences. Fi and ifg depend 
on an arbitrary factorization scale fi in such a way that the physical decay rate is 
independent of fi. In order to avoid large logarithms of m-f/^ in Fi, /i of 0(mx) 
should be used. The factorization theorem for P-wave charmonium decay contains a 
nonperturbative parameter H^ analogous to ifg. However a production process must 
be used to determine H'^ phenomenologically f^. 

In the usual nonrelativistic quark model. Hi can be expressed in terms of the 
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P-wave color-singlet cc wave function: 



H, ^ A . 15 MeV, (2) 

ZTT 



c 



where the numerical estimate was obtained in Ref. [0 from measured decay rates of 
the Xci and Xc2- 

ifg cannot be rigorously expressed perturbatively in terms of Rp, since it accounts 
for radiation of a soft gluon by a color-octet cc pair. The scale dependence of -ffg(/i) 
is determined by the following renormalization-group equation (to leading order in 

d 16 
^dfl^'^^^'' (3) 



which is readily integrated. For example 



16 , / asjfJ'o) \ 16 faJjricA 



Hi (4) 



27/33 V"s("^c)/ 27/34 V"^("^b)/. 
(for /io < mc), where /3„ = (33 — 2n)/6. If H'^^hq) is neglected in the limit of large 
one obtains Hg{mh) ^ 3 MeV, using asifJ^o) ~ 1 While one might not expect the 
physical value of rrij, to be large enough to neglect H^{^q), an estimate for H^{mb) 
obtained in Ref. ||^ from experimental data on B meson decays is consistent with 
the above result. 

A calculation of f i and f § in Eq. (|l|) each to leading order in as can be obtained 
from a calculation of the infrared divergent width r^iv for T — * cc('^P7) + ggg, where 
the cc pair is in a color-singlet P-wave state: 



rdiv(T cc{^Pj) + ggg; /io) = 
20a^^ Gj 



Hi. (5) 



J-'ij and jFg are dimensionless infrared-finite form factors, /io is an infrared cutoff on 
the energy of soft gluons, and fi is an arbitrary factorization scale [the /i dependence 
of J^ij exactly cancels that of the explicit logarithm in Eq. (^]. 

The constants in Eq. (|^) include a color-factor of 5/216 and phase space fac- 
tors, including 1/3 for T spin-averaging, and 1/3! for the phase space of the three 



indistinguishable gluons [cf. Eq. (|1^) below]. Gj is related to the usual S'-wave bb 
nonrelativistic wave function: 

Gj ^ ^ 108 MeV, (6) 

where the numerical value is obtained from the electronic decay rate of the T 0. 

The hard subprocess rates of Eq. (|ID are identified from Fdiv by using the per- 
turbative expression for the infrared divergence in H^, obtained from Eq. (|^) by 
neglecting the running of the coupling 

^8(/^)~T^«.lnf-Wi. (7) 
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Thus: 



and 



fi(T ^ ccCPj) + ggg; /i) = — J-^X/i), 



f8(T ^ cc{'S,) + gg) = t^^^s- (9) 

Note that f i is suppressed by 0{as) compared to f §. However, the nonperturbative 
parameters Hi and which accompany these subprocess rates in Eq. (|1]) are inde- 
pendent, hence agHi need not be small compared to H'^ [0, We therefore proceed 
to calculate f i and f § each to leading order; all further corrections to P-wave pro- 
duction are then guaranteed to be suppressed by at least one power of compared 
to what is included here. 

In order to extract J^u and JFg individually, it is necessary to explicitly identify 
the infrared logarithm in the calculation of Fdiv This can be done analytically, as 
described in the following. 

There are 36 0{al) diagrams contributing to Fdiv One of these is shown in 
Fig. |I[ All infrared divergences are associated with a gluon that is radiated from a 
charm quark line; in Fig. |I| this gluon carries four momentum ki and polarization ei. 
Define the invariant amplitude A4j{2, 3; 1) corresponding to the sum of all Feynman 
diagrams where gluon "1" is radiated from the charm quark line. The amplitude 



is readily computed using expressions for S- and P-wave QQ currents given in Ref. 

M(2'^--\) = mrm^B,{2,3)C%l) 

' ' '~ [ik2 + k,) ■ k,][{k, + k,) ■ k^][{k^ + ^3) ■ k,] kl {k ■ k,r ^ ^ 

where 

e^S^(2, 3) = {e4 ■ 62 [-^4 ■ he^i ■ k2eo ■ /C4 - ^2 ■ ^sCs ■ ^460 ■ ^2 - ^4 ■ ^3^2 ■ k^eo ■ €3] 
+ eo ■ £3 [^4 ■ ^3^4 ■ k2e2 ■ k3 + k2- k3e2 ■ k^e^ -k^-k^- k2ei ■ k3e2 ■ k^]} 
+ {2 ^3} + {3 ^4} (11) 

(eo is the polarization of the T), and 

e4^C^^o(l) = [d ■ £4/^1 -k^-ei- k^e^ ■ ki] (^mj + k ■ k^ - kj^ , 

e4;,Cj=i(l) = ^m^kleap^se°e^elk(, (12) 

e4/.C^=2(l) = (kl ■ A;4e?ef + k^kle^ ■ - fc^efei ■ k^ - k^e^e^ ■ ki) e"^. 

is a spin-1 polarization vector and e"^ is a spin-2 polarization tensor. 

For convenience the virtual gluon is labeled in Eqs. (|TD|)-(|I2]) by polarization €4 
and momentum ^4 (^4 = P — A;2~^3 = ^ + ^i)- Terms which vanish due to the on- 
shell conditions ■ /cj = (i = 1, 2, 3) and eo ■ P = have been dropped. The charm 
quark current Cj(l) is symmetric under interchange of labels 1 and 4 (up to terms 
which vanish due to the on-shell conditions). The bottom quark current P^(2,3) is 
explicitly symmetric under interchange of labels 2, 3 and 4. 

The overall factors in Eq. (|^) are such that: 

3 / <i|4.4l E [-^5(2. 3; 1) + 2MA2, 3; l)Mj(l, 3; 2)] , (13) 

spins 



^Overall factors in the quark currents including couplings, color amplitudes, and wave functions 
have been accounted for in Eq. 
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where denotes (infrared-cutoff) n-body phase space, normahzed according to 



i=l 



The factor of 3 on the right hand side of Eq. ([T3| ) accounts for symmetrization of 
A^j(2,3; 1) under gluon label interchanges 1^2 and 1 3, taking account of the 
symmetry in the three gluon phase space. 

The infrared divergence comes entirely from the first term in square brackets in 
Eq. (p!3|), and is due to the P-wave charm quark propagator Ijik.kiY in Eq. (plQl). It 
is therefore advantageous to organize the four-body phase space integral in Eq. (|13D 
by taking the invariant mass of the Xcj and gluon "1" as one integration variable 

d[$4] = / d{kl,) / d{kl^) 

JO Jmj+2^o™x 

X $2[P ^ k23, /cij '^2[k23 h, h] ^^21^^ h, k], (15) 

where 77123 = ^23- Note the infrared cutoff /iq on the energy of gluon "1" in the rest 
frame of the Xcj- 



The infrared logarithm on the right-hand side of Eq. (|T^) can now be identified 
analytically by observing that P^(2,3)Cj(l) in Eq. ( PUj) is given by a sum of terms 
each containing exactly one factor of /ci, if ^4 = P — A;2 — ^3 is used to eliminate the 
virtual gluon momentum. With this convention, one has 

^A,5(2.3;l) = 2^(^t^|^. (16) 

spins V 1 z 

where ki appears explicitly in the function 7j(fci; P, k, k2, k^) only in the combination 
ki/k ■ ki. 

J-'s is then given in terms of a manifestly infrared-finite three-body phase space 
integral, taking account of the fact that $2(^ix ~^ ki, k) = \k ■ ki/kf^J dfll^, where 
is the center-of-mass solid angle of the two body system: 

27tt /■(™T~'Ttx)^ 

(2J+l)J-8 = d{kl,)<i>2[P^k23,k] 

X 

X $2[/c23 ^ k2,k3] J dni^-fj{ki;P,k,k2,k3), (17) 



where 

ki= lim (18) 
k-ki^o k- ki 

The finite four- vector ki is readily expressed directly in terms of and fi^^. An 
expression for J^u can be obtained from Eqs. ( p!3D and (|l^) by analogy with the 
identity / dxf{x)/x = /(0)lnx + / dx[f{x) — f{0)]/x. 



The contraction of currents and sum over polarizations in Eqs. (0) and (|13|) were 
performed symbolically using REDUCE H (leading to lengthy expressions, particu- 
larly for J = 2). The Xcj spin sums were done using (see e.g. Ref. 0): 

k k 

(^[i^v — Qfiu -r 2 ~ ' 

e/.i.e„/3 = I [P^aVup + V^pVua] - IV^uVap- (19) 

e 

The phase space integrals were evaluated numerically using VEGAS [[1^]; modest 
integration grids are found to give very good convergence. 

The fact that jFg should be independent of J provides a stringent check of these 
calculations, given that the three currents Cj have very different structures [cf. Eq. 
(|T2|)]. This was verified explicitly in numerical calculations of Eq. (p!7D, to better 
than a few tenths of a percent for all m^/rrir on a modest integration grid. 

Figure |^ shows the numerical results for jFg over a range of hypothetical meson 
masses. In Fig. |^ results for J^ij{fi) are shown using a factorization scale /i = m^. 

The available experimental data on charmonium production in T decay is for the 
J/^: 

= (1.1 ± 0.4) X 10-3 CLEO [11], 

< 1.7 X 10-3 Crystal Ball pl, (20) 

< 0.68 X 10-3 ARGUS [lU 



5exp(T^ J/V^ + X) 



An upper bound on H'^ can be extracted from this data by computing the "indirect" 
production of J/ip due to the Xcj states. Assuming that radiative cascades from Xci 
and Xc2 dominate, with branching fractions -Bexp(Xci iJ/i^) ~ 27% and -Bexp(Xc2 ^ 
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'^J/ip) ~ 13% 0, the results presented here give: 

H^(,„,)«{lk£(I^||±£;^i^,,4)MeV. ,21) 

The first number in brackets above comes from the color-octet subprocess rate Fg, 
and the second number from the color-singlet rate Fi. The experimental value for 
the total width FtotC^) ~ 52 keV was used, along with asirny) ~ 0.179 [0], and 
the values of Hi and Gj given in Eqs. (g) and (|]). 

Equation yields the bound H'^{mi)^2'b MeV using the ARGUS upper limit, 
which is consistent with the other measurements. This bound is considerably larger 
than an estimate H^{mh) ^ 3 MeV based on B meson decays IQ,*^ although a 
calculation of next-to-leading order QCD corrections to the color-singlet component 
of i? — > Xcj + ^ is required before an accurate determination of ifg can be made in 
that case |0. 

This raises the possibility of significant direct production of J/ip in the decay of 
the T, unless the branching fraction turns out to be considerably smaller than the 
ARGUS bound. Mechanisms for direct T ^ J/ip + X in perturbative QCD were 



first discussed in Refs. ||14[ and ||T5l. The direct production rate is suppressed by 
0(a;^) compared to the P-wave color-octet production mechanism considered here. 
However, the nonperturbative matrix elements which enter into P-wave production 
are of 0{v^) relative to the corresponding parameter for 5- wave production, where v 
is a typical relative velocity of the quarks. Moreover, there are many channels which 
contribute to direct production. 

The full 0{a^) perturbative QCD amplitude for direct T J/i/j + X was recently 
evaluated in Ref. [16[, corresponding to one- loop diagrams ioT T ^ J/ip + gg, and 



tree diagrams foT T J/ip + gggg- The 0{a1a'^) electromagnetic amplitude for the 
two gluon decay mode was also evaluated. Unfortunately, only a crude estimate of the 
required phase space integrations was made in Ref. [1^ (there is a costly convolution 
^From Eq. Hg{ii) increases by only « 0.3 MeV in the evolution from /i = mb to i^i = mx- 
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with a numerical calculation of the loop integrals for T ^ J/ip + gg). Nevertheless, 



the calculation of Ref. [0 suggests a branching fraction for direct production of a 
few X 10"''. This would lead to a considerable reduction in the bound on ifg extracted 
from Eqs. and (|2l|). 

To summarize, a complete calculation was made of the leading order rates for 
^ ~^ Xcj + X, through both color-singlet P-wave and color-octet S-wave cc subpro- 
cesses. Experimental data on J/ ip production was used to obtain an upper bound on 
the nonperturbative parameter ifg, related to the probability for color-octet S-wave cc 
hadronization into P-wave charmonium. This work adds to the limited information so 
far available on the color-octet mechanism for P-wave quarkonium decay and produc- 
tion 1^, ^. These investigations provide new information on a nonperturbative sector 
of QCD that has largely been neglected in previous studies of heavy quarkonium. A 
quantitative estimate of H'^ is phenomenologically important since this parameter is 
required as input for the calculation of a variety of processes. Improved experimental 
data, and a definitive calculation of the direct J / ip production rate along the lines of 
Ref. [1^ , would allow for an accurate determination of from the results presented 
here. 
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Figure 1: One of the 36 0{al) diagrams contributing to T — > cc{^Pj) + ggg. 



Figure 2: Color-octet form factor as a function of m^/my. 



Figure 3: Color-singlet form factors J^u as functions of m^/mr'- J = (short-dashed 
line), J — 1 (long-dashed line), J = 2 (solid line). The form-factors were evaluated 
using a factorization scale /i — mr- 
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